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Mai or ana- fermionic quantum computation (MFQC) was proposed by Bravyi and Kitaev (See 
Ref.[l|), in which a fault-torrent (non-topological) quantum computer built from Majorana fermions 
may be more efficient than that built from distinguishable two- state systems. However, till now 
people don't know how to realize a MFQC in a physical system. In this paper we proposed a 
possible realization of MFQC. We find that the end of a line-defect of p-wave superconductor or 
superfluid on a honeycomb lattice will trap a Majorana zero mode, which will become the starting 
point of MFQC. Then we show how to manipulate Majorana fermions to do universal MFQC, which 
possesses unique possibilities for high-level local controllability, individual addressing, and readout 
of the quantum states of individual constituent elements by using timely cold-atom technology. 



Introduction. — Quantum computer makes use of the 
principles of quantum physics to enhance the compu- 
tational power beyond what is attainable by a tradi- 
tional computer. Various designs have been proposed to 
build a quantum computer, such as manipulating elec- 
trons in a quantum dot or phonon in ion traps, cav- 
ity quantum electrodynamics, nuclear spin by nuclear 
magnetic resonance techniques. However, the realiza- 
tion of low decoherence condition in these schemes is a 
great challenge to achieve a true quantum computer. Be- 
cause the noise and decoherence are presumably caused 
by local interactions, people try to avoid them by en- 
coding quantum information non-locally. This leads to 
so called topologically protected qubit[2|-[3. In partic- 
ular, the topological quantum computation is proposed 
based on the manipulation of non-Abelian anyons^5ji9j]. 
The possible candidates may be the fractional quantum 
Hall state at filling factor = 5/2 in ultra high- mobility 
samples plol [Tlj , the two dimensional (2D) chiral Px + 
ipy superconductors [HI, and the s-wave-superconductor- 
topological-insulator systems [l3|, Yet, people cannot 
realize universal quantum computation just by braiding 
Ising-type non-Abelian anyons. In addition, among these 
approaches, accurate manipulation of single anyon re- 
mains a major difficulty and new techniques are expected 
to overcome this drawback. 

In this paper we turn to an alternative type of 
fault-torrent quantum computation - Majorana- fermionic 
quantum computation (MFQC). Major ana- fermionic 
quantum computation was proposed by Bravyi and 
Kitaev [H, in which the qubits are characterized by 
fermion's parity. A quantum computer built from Majo- 
rana fermions may be more efficient than that built from 
distinguishable two-state systems. Although MFQC is 
not the topological quantum computation, the quan- 
tum information with Majorana fermions also intrinsi- 
cally immunes to decoherence. Thus it is possible to 
do universal quantum computation by manipulating the 
two-Majorana-fermion interaction (7172) and the four- 
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Majorana- fermion interaction (71727374). We find that 
the end of a line-defect of p-wave superconductor (SC) or 
superfluid (SF) on a honeycomb lattice will trap a Majo- 
rana zero mode, which will become the the starting point 
of MFQCp^]. Then we show how to manipulate Majo- 
rana fermions to do universal MFQC and also discuss the 
errors. In the end, we design an experiment scheme for 
the realization of the p-wave SF in a Bose- Fermi mixture 
on an optical honeycomb lattice. 

P-wave superconductor / superfluid for spinless 
fermions on honeycomb lattice: We start from a 
spinless fermion model on a honeycomb lattice with 
nearest attractive interaction as 



H = —t^^^clcj—t^ c\cj^h.c. — pL^y^^fij — U ^y^^hjhj, 

(1) 

where t {t') and U denote the strength of nearest (next 
nearest) neighbor hopping and the nearest neighbor at- 
tractive interaction respectively, fi is chemical potential, 
and hi = cJq denotes the particle number operator. 

As increasing the interaction strength, the system 
turns into a p-wave SC/SF ground state in the region 
of 3.36t < U < 1.12t for the case of = at zero 
temperature • In the limit oi t' <C t, the p-wave 
SF order is still stable. For the p-wave SC/SF, the 
p-wave pairing order parameter is given as Aj^j+ai = 

c]c^^^^^ = 0, where denotes a vector that connects 
nearest neighbor i site and i + Bq, site as ai = |(1, >/3), 
a2 f(l,-\/3), as a (-1,0), bi ai - as, bi 
^2 — as. The lattice constant is defined as a = 1 in this 
paper. In Figdja), along red links, the SC/SF order 
parameter is finite; while along black links, the SC/SF 
order parameter is zero. The effective Hamiltonian then 
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FIG. 1: (a) The illustration of honeycomb lattice of two in- 
terpenetrating triangular lattice. A and B denote A and B 
sublattice. The red links represent the paring bonds, on which 
the pairing order parameter is finite; (b) The order parameter 
A vs. the interaction strength U for the case of t' = O.Olt at 
half-filling; (c) The order parameter A vs. the temperature 
T for the case of = O.Olt, U = 5.24t at half-filling. 

becomes 

{ij) {{ij)) jeB 

jeB i i 

(2) 

At zero temperature, there exists a first order phase 
transition between semi- metal and the SC/SF phase, 
where Uc — 3.23t when f = O.Olt at half filling as 
shown in Fig. 1(b). In weakly interaction region, U <Uc-, 
the ground state is a semi-metal. When the interaction 
strength is larger than Uc-, in the region of 3.23t < U < 
7.12t, the ground state turns into a p-wave SC/SF state. 
In the following parts, we focus on the case ofU = 5.24t, 
of which the SC/SF order parameter A 1.34t at zero 
temperature. As the temperature increases, the order pa- 
rameter decreases until it disappears at the critical tem- 
perature T* = l.lOt that corresponds to breaking the 
Cooper pairing. The order parameter A versus temper- 
ature T is plotted in Figdjc) for the case of t' = O.Olt, 
U = 5.24t at half-filling. 

For the p-wave SC/SF on honeycomb lattice, there ex- 
ists an energy gap for the Bogliubov quasi-particles. For 
example, for the case of U = 5.24t, the energy gap AEg 
is about 2.0t. Remember that for the p-wave (p^^-wave 
or py-waxe) SC/SF state on a square lattice, the spec- 
trum of quasi-particles is always gapless. And it is ob- 
vious that this p-wave SC/SF on honeycomb lattice is 
not a "topologically ordered" SC/SF. However, there ex- 
ist Majorana edge states. Now we consider the effective 




FIG. 2: (Color onhne) Edge states of p-wave SC/SF on hon- 
eycomb lattice for the system with armchair open bound- 
ary. The parameters are chosen as A = 1.34t, t' = O.Olt, 
/J — 0.0136t. Along red links, the SC/SF order parameter is 
finite; while along black links, the SC/SF order parameter is 
zero. 



Hamiltonian in Eq.(|2]) under periodic boundary condition 
along y-direction and opened boundary condition (arm- 
chair boundary) along x-direction (the axis is defined in 
Figl2fa)). Now the p-wave SC/SF on a honeycomb lat- 
tice behaves like coupled ID p-wave SC/SFs along x- 
direction, where the dangling Majorana fermions exist at 
the boundary [igI]. Fig l2] shows the numerical results of 
the edge states for the case of A = 1.34t, t' = O.Olt, 
/i = 0.0136t. From FiglSl one can see that the Majorana 
edge states always have a flat-band without hybridizing 
with each other. 

Majorana zero modes around vacancies: Now we study 
•the properties of a lattice defect for the p-wave SC/SF 
on honeycomb lattice in the presence of an on-site po- 
tential at site io, ^eff ^eff + "^io^o^^o- After solving 
the Bogolubov-de Gennes (BdG) equations with a single 
lattice defect on a 50 x 50 lattice with periodic boundary 
condition, we obtain two localized wave-function with en- 
ergy ±e around the lattice defect. In the unitary limit, 
the lattice defect becomes a vacancy shown in Figl3l^a), of 
which we have an infinite on-site potential, i.e., Vi^ ^ oo. 
For the system with the particle- hole symmetry, /i = 0, 
t' = 0, two localized states turn into a pair of Majorana 
modes with exact zero energy [isj, one Majorana mode 
at the left side of the vacancy, and the other at the right 
side. See the results in Figl3fa): there are two Majorana 
modes, 7^, 75. Away from the unitary limit, we have a 
finite on-site potential. As shown in FigjSl^b), the energy 
levels of the two localized modes versus potential energy 
is plotted. From Figl3fb), we find that as Vi^ ^00, 
the energy of the localized state turns into zero. 

Next, we consider the case of a line of lattice-defect 
with multi- vacancy along the SC/SF paring direction as 
shown in Fig|4](a). Now the Hamiltonian becomes i^eff 

^eff + Xl^i "^ifccjfc^^fc- -^y numerical calculations on a 
50 X 50 lattice, from FigHJ^a), we find that there also 
exist two Majorana modes localized at both ends of the 
line-defect, 7^, 75. Now, the end of a line-defect can 
be considered as the boundary of ID p-wave SF. Thus 
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FIG. 3: (a) The lattice vacancy is represented by shaded- 
dashed-circle, and two green dots denote the two Majorana 
modes 7a, 76- Along red links, the SC/SF order parameter 
is finite; while along black links, the SC/SF order parameter 
is zero, (b) Half of the energy splitting Siji versus on-site 
potential Vi^lt. The detailed explicit Hamiltonian form is 
shown in Eq.([2|)- Related parameters are chosen as A = 1.34t, 
= 0, /i = at zero temperature. 
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FIG. 4: (a) The illustration of a line-defect. Two Majorana 
modes 7a, 7b are denoted by green dots are localized at ends 
of the line-defect, (b) The illustration of line-defect created 
by an addressing line-laser in optical lattice. Along red links, 
the SC/SF order parameter is finite; while along black links, 
the SC/SF order parameter is zero. 

each end of the line-defect traps a dangling Majorana 
zero mode. When the line-defect is long enough, the two 
Majorana modes decouple and have zero energy. 

Universal Majorana Fermionic Quantum Computa- 
tion: Then we consider the couplings between Majorana 
zero modes of vacancy-lines. For a general case, the uni- 
versal Hamiltonian of Majorana zero modes is given by 

i^Majoaran = i ^ Jijialib + i ^ Jiflihlja (3) 
^ ^ Jij ^ lia^Jj alibi jbi 

where Ji denotes the coupling between two Majorana 
modes of a line-defect i, J-^J^ denotes the coupling be- 
tween two Majorana modes of two different line-defects 
i, j, and J^^^ denotes the strength of the four-Major ana 



particle interaction of two different line-defects i, j. 7^^ 
denotes a Majorana fermion at the left side of a line- 
defect i and 7^5 denotes a Majorana fermion at the right 
side of a line-defect i. 

Therefore, the two Majorana modes (7^^, Jib) of a line- 
defect form a physical fermion's parity qubit, |0) • (quan- 
tum state with even fermion's parity) and |0)^ = 
(quantum state with odd fermion's parity), respectively. 
We label a pair of Majorana edge states (7^^, Jib) by a 
complex fermion as jia = (ci + cj^ , 7^5 = -i (^Ci - 4) • 
By representing the Majorana fermions into complex 
fermions, the universal Hamiltonian in Eq.(j3j) of coupled 
Majorana zero modes becomes 

i^Majoaran = ^Ji (^4^^ " l) (4) 

{ij) 

-Y.j^ H^^-i) H^.-i)- 

We firstly calculate J^, the coupHng between two Ma- 
jorana fermions at the ends of a single lattice defect (a 
shortest line-defect with a defect). In the unitary limit, 
we have two degenerate zero modes (7^^, 7i6) around the 
lattice defect for the SC/SF with particle- hole symme- 
try (/i = 0, = 0). Away from the unitary limit of the 
SC/SF, Vq 7^ 00, we have the finite energy splitting of the 
Majorana modes as Si 7^ 0. Then we can identify half of 
the energy splitting £i/2 to be the coupling between two 
Majorana modes around this lattice defect as Ji =6^/2. 
When the energy splitting > (< 0), the system en- 
ergetically favors the empty (occupied) state. We also 
calculate the energy splitting of two Majorana fermions 
of a line-defect with multi-defect (the number of the de- 
fect is more than two) and found that the energy splitting 
of the two Majorana modes is always zero. In this sense 
from Figini^b), to tune the coupling Ji (that is £i/2) of 
a line-defect between two Majorana modes (7^^, lib)-, we 
can first shrink the line-defect into a single vacancy and 
then change the on-site potential Vq from infinite to a 
finite value. 

Next we calculate JI- , the coupling between two Ma- 
jorana modes at the ends of two different line-defects. 
Now we consider two line-defects as shown in Fig. 5 (a). 
When two line-defects are well separated, we have four 
Majorana modes 7^^, 7i6, Tja, Ijb with exact zero energy. 
Now we move the two line-defects nearby. The wave- 
functions of the two Majorana zero modes 7^5, 7^^ will 
overlap and give rise to a small energy splitting /S.Eij 
from exact zero energy level. The tunneling splitting 
/S.Eij is about J^^^ ex e~^^9^*'j where Li^j is the dis- 
tance between right-site of z-th line-defect and left-site 
of the j-th line-defect. In Fig|5l^b), we plot the results 

of J^^^ versus Li^j for the case of t' = O.Olt, fi = 0.013t. 
For this case Li^j » a, we get zero tunneling ampli- 



FIG. 5: (a) The illustration of two line-defects. The two 
Majorana fermions 7ia, Jib are denoted by the green dots 
that localize at ends of the line-defect 1 and the line-defect 
2, respectively, (b) The energy splitting AE/t versus the 
distance between the ends of two line-defects as shown in (a). 
The parameters are chosen as A = 1.34t, t' = O.Olt, /i = 
0.013t. Along red links, the SC/SF order parameter is finite; 
while along black links, the SC/SF order parameter is zero. 



tude as Jij ^ 0. In this sense, to tune the coupling J-j ^ 
(that is AEij) of two line-defects between two Majorana 
modes (7^5, jja)^ we can just change the distance of two 
line-defects. 

Thirdly we calculate the strength of four-Majorana 
particle interaction jj^j^ between the four Majorana 
modes jia^ Jih-, Jja^ Ijh around two line-defects. Here 
we also consider two line-defects in the unitary limit. 
When two line-defects are well separated, there is no 
four-Majorana particle interaction at all. When two line- 
defects are nearby, the four-Majorana particle interaction 
could be induced by the nearest neighbor attractive inter- 
action along un-pairing direction as —U^-hifii-^as- As 
shown by FiglHl there are four Majorana fermions around 
the two line-defects 7^^, 7i6, 7i+a3,a, 7i+a3,6. The two zero 
modes 7^^, Jib form a complex fermion q by the defini- 
tion jia = (ci + 4) ' 7i6 = -i (ci - 4^ and the other 
two zero modes 7^+a3,a, 7i+a3,6 form another complex 
fermion cj by the definition 7i+a3,a = (Q+a3 +4+a3) ^ 

7i+a3,6 = -i (^Q+a3 - 4+a3)- Then after considering the 
nearest neighbor attractive interaction —U^-fiifiij^s,^^ 
we will get residue interaction term between the two com- 
plex fermions (q, Q+a3)- For this case, the strength of 
the four-Majorana particle interaction term is given by 



-J. 



(4) 

i,i+a3 



i 



U 



(5) 



In this sense, to tune the four-Majorana particle inter- 
action J^"^^ of two Hne-defects between four Majorana 
modes {jia-, Jib-, Ija-, Ijb)^ we can just move two par- 
allel line-defects nearby. However, for this case, the 
two-Major ana couplings between (7^^, 7i+a3,a) and (7^5, 



FIG. 6: The illustration of four Majorana modes induced by 
two defects that can combine into two complex fermions (q, 
Ci+a3)- Along red links, the SC/SF order parameter is finite; 
while along black links, the SC/SF order parameter is zero. 



7i+a3,6) are still finite, Jljj^^^ ^t^O. 

By freely tuning all terms in Eq.(|l]), we may make 
the following two universal gate set in fermionic quan- 
tum computation (up to phase factor) as e^^^^j and 
g*f 71727374 p. design e^^'""^^^ gate we can shrink 

the line-defect into a single vacancy and then fix the 
on-site potential Vq under a given time-period as At = 
8J-^(Vo) ' design e^^^^^^"^ gate we can control the 

given two line-defects nearby under fixed time-period as 
= —(^ ; 3) To design e^^^^-^^^^^+-3, -7^+^3,5 gate 

we can control the given two line-defects nearby under 
fixed time-period as At = — — [22|. Thus, in princi- 

pie, people are capable of doing arbitrary unitary trans- 
formation on the protected subspace by controlling the 
vacancies. 

Errors: For the p-wave SF on honeycomb lattice, one 
always has a large energy gap AEg for the fermions and 
a very tiny energy splitting AE of the nearly-degenerate 
Majorana modes, i.e., AE <C ^Eg. Based on this con- 
dition {AE <^ we may ignore the excited states 
with E > AEg and consider only the Majorana modes as 
the protected sub-space. In this sense, we must consider 
the case at fairly low temperature ksT <C AEg. Now 
there exist few excited quasi-particles, n/ ~ ^-AEg/ksT ^ 
In addition, the random local potentials will lead to the 
energy splitting of the Majorana modes. Thus when we 
move the vacancies nearby to do unitary operations, the 
random energies of the Majorana modes will lead to er- 
rors. By numerical calculations, we find that the energy 
splitting by random potential is very tiny. So we may 
ignore this type of errors. 

To store quantum information by the Majorana zero 
modes of the line-defect, we must use a long line-defect 
to separate the Majorana modes. Thus, the distance be- 
tween two Majorana modes are set to be much bigger 
than lattice distance a as min Li j » a. For this case, 

we get zero tunneling amplitude as J!>j {Lij) 0. Now 
the system truly immunes to decoherence. 

Possible realization in cold atoms: Now we propose 
an experiment scheme for the realization of the spin- 
less fermion model on a honeycomb lattice with near- 
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est attractive interaction given in Eq.([T]). To realize this 
model, a Bose-Fermi mixture on optical honeycomb lat- 
tice of cold atoms may be a possible candidate. As 
given in Ref.|23|], the bosons with unit filling form a 
Mott insulator and the spinless fermions with half fill- 
ing form a semi-metal. When the on-site repulsive inter- 
action between bosons and fermions Ubf is larger than 
boson-boson repulsion Ubb i-e., Ubf > Uhb-, a spinless 
fermion and a bosonic hole bind into a composite spin- 
less fermion created by = ph^ with /''', h being the 
fermionic creation and the bosonic annihilation oper- 
ators, respectively. The effective nearest attraction is 
given by U/t = 2 {Ubf /U^h — ^) ^ and the effective near- 
est (next nearest) hopping t = 2tlj/Ubf {f = 2t'^^/Ubf)^ 
where tbf (tbf) denotes the tunneling strength for the 
bosons and fermions of the hopping between the nearest 
(next nearest neighbor) neighbor lattices (here we have 
supposed that tb = tf = hf^ t'^ ~ ^'f ~ ^bf 
plicity). See detailed calculations in the supplementary 
materials. In experiment, for example, ^^^Yb— ^^^Yb 
or ^Li— ^Li mixture, the nearest neighbor effective at- 
tracting interaction U may be tuned employing well- 
characterized homonuclear and heteronuclear Feshbach 
resonances (23 - [26| . Thus in cold atom system, people 
have ability to tune the parameters in Eq.(l) including 
f/t, U/t. 

On the other hand, recently, in cold atoms people 
have made great progress on the manipulation of a sin- 
gle atom on an optical lattice p7l-[29|. A tightly focused 



dimple laser beam with an accuracy about 0.1a can be 
applied onto individual lattice sites to generate a vacancy 
in the optical lattice or line-vacancv |27l - [29j . By moving 
the addressing laser beam as shown in Fig. 4(b) from one 
site to another, people are able to accurately manipulate 
the position and the length of a line-defect on a two- 
dimensional optical lattice. Therefore two dimensional 
ultrocold atomic Fermi gases provide an ideal platform 
for creating and manipulating Majorana fermions. 

Conclusion: We have proposed a feasible scheme to- 
wards designing a quantum computer, which incorpo- 
rates intrinsic fault tolerance by controlling the Majo- 
rana fermions accurately via controlling addressing laser 
beams on optical lattice. In particular, we find that the 
lattice defects here plays the role of a Majorana fermion. 
Therefore, by creating and manipulating lattice defects 
accurately via addressing laser beams, we can do quan- 
tum computation fully optically. So we call this approach 
Majorana- fermionic quantum computation'''' to empha- 
size the difference between the topological quantum com- 
putation by braiding non-Abelian anyons. 
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